Universal high-frequency transport in perfect photonic crystals 
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The light scattering in the periodic dielectric cylinder array is studied. We analytically calculate 
the diffusive-ballistic transport crossover and find the weak localization superimposing on it. Possible 
experimental observations are analyzed. 
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Introduction — Experimental studies [H, 0] are offering 
new boost to interest in transport of optical systems. A 
traditional subject is the light scattering in fully ^] or 
partially [1, 3] disordered media. Irrespective of the ab- 
sence/presence of the periodic background disorders turn 
out to be essential in both kinds of materials. There, the 
Anderson localization of light originates at the multiple 
disorder scattering. The latter leads to the enhanced 
constructive interference between two counterpropagat- 
ing optical paths and, thus, suppresses the light diffusion. 

For perfect periodic dielectric materials (namely pho- 
tonic crystals) the low-frequency transport is protocoled 
by sample-specific band structures, the analysis of which 
is nowadays a well-established industry. By contrast, the 
high-frequency transport in perfect photonic crystals is a 
largely unexplored subject which might find practical ap- 
plications. There, localization and periodicity may inter- 
play with each other. Theoretically, no localization sig- 
natures have been found for periodic quantum mapping 
which might suggest that localization is incompati- 
ble with perfect periodicity. Contrary to this, for realis- 
tic quantum particle motion in some perfect crystals, the 
weak localization correction to the classical diffusion coef- 
ficient is analytically shown Q . Experimentally, to detect 
localization signals in periodic structures, if any, turns 
out to be extremely difficult. Indeed, in usual systems 
(e.g. electrons) ubiquitous particle interactions compli- 
cate analysis of wave interference effects such as (weak) 
localization. Although this problem is circumvented in 
optical systems, traditional transmission/reffection mea- 
surements do not allow directly probing the dynamical 
(bulk) diffusion coefficient. 

The recent experiment [1] offers the very opportu- 
nity to study the diffusive-ballistic transport crossover. 
The experimental scope consists of two crucial steps: (i) 
to fabricate the dielectric material which is uniform in 
the longitudinal direction (but might be arbitrary in the 
transverse plane) and, (ii) to launch a probe laser beam, 
parallel to the longitudinal direction, into the material 
(Fig. [T]). Then, the dynamical process of the dispersion 
in the transverse plane is recorded by the intensity profile 
at different propagating distances, which allows the di- 



FIG. 1: A narrow input probe beam (green spot)- 
perpendicular to the plane-is launched into the dielectric 
cylinder array and expands. The "time" axis points inwards. 



rect measurement of the dynamical diffusion coefficient. 
It is the purpose of this letter to analytically study the 
crossover in perfect photonic crystals and to show that 
weak localization superimposes on this crossover. Fur- 
thermore, we examine the possibility of confirming these 
predictions within this experimental scope. 

Qualitative results based on band structure — To be spe- 
cific we will consider the following photonic crystals. Di- 
electric cylinders of a radius , uniform in the "longitu- 
dinal" (z-) direction, are periodically embedded into the 
host material of a refractive index Uq forming a triangu- 
lar lattice in the "transverse" {x-y) plane with a lattice 
constant a (Fig. [T]). For this lattice let us choose the 
coordinate system as {x,y) — 1 + r , where 1 G 05 (the 
Bravis lattice) is chosen to be the cylinder center, and r 
specifies the position inside the unit cell. A probe beam 
of frequency fio , parallel to the z-axis, is launched into 
the host dielectric material, then expands as propagates 
along the z-direction. The transverse wave number k± is 
inversely proportional to the initial width of the beam, 
which is narrow enough such that kj^ <^ min(a— 2ro, tq) . 

Then, the propagation of the scalar field (for sim- 
plicity) denoted as E{x, y, z) inside the photonic crys- 
tal is described by the Helmholtz equation (The light 
velocity in the vacuum is set to be unity.): [V^ + 
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V,Qn'^{x, y)]E{x, y,z) = with n{x, y) the refractive index 
profile. Introduce the slow-varying envelope il^i^i Hi z) 
such that E{x, y, z) = ^{x, y; z) e^^°^°^ and insert it into 
the Helmholtz equation. In the paraxial limit \d'^ip\ ^ 
f^o'^olc^zV'l the second order derivative may be ignored 
[H, consequently we arrive at the Schrodinger equation: 



{dl + dl)-V{x,y)\^{x,y;z) = 0, 



where V{x, y) 



(1) 

This shows that 



the beam propagates in the manner fully analogous to the 
two-dimensional quantum particle motion in the "poten- 
tial" V{x,y) , and z plays the role of the usual time. 

In periodic potentials exists the band structure e„(k) , 
where n is the band index and k is the quasimomcntum 
in the first BriouUin zone (B.Z.). As to be detailed below 
(for a large class of potentials 6]) the ray optics inside 
a unit cell is ergodic characterized by the positive Lya- 
punov exponent A . Following the well-known Bohigas- 
Giannoni-Schmit (BGS) conjecture, then, given k (away 
from symmetry points or lines) the high-frequency bands 
display universal level fluctuations ■ 

If a single band is excited the motion of (Bloch) pho- 
tons is ballistic. However, because the beam initially 
is narrow enough it excites a number of high-frequency 
bands centered ai e — fc^/(2r2o"-o) with an amount of 
^ A/ A ^ 1 , where A is the mean level spacing. At early 
times the rays are successively bent by dielectric cylinders 
which, equivalently, may be viewed as that Bloch pho- 
tons experience considerable inelastic scattering. Con- 
sequently, the beam dispersion (expansion), defined as 
a^{z) = J dxdy{x^ + j/^)|?/'(a;, y; z)p , is dominated by 
the interband transition factor: e*'''^"^''-'"'^"'^'^^^^ . Then, 
the large number of excited bands allows us to average 
this factor with respect to the universal level correlator, 
which is ~ {en(k) — e„'(k)}~^ fs', Consequently, the 
normal diffusion, i.e., cr'^{z) ^ z is recovered for early 
times A^^ < z <C A^^ . At later times z > A^^ indi- 
vidual bands are resolved and the interband transition 
becomes negligible. The transport, in turn, becomes bal- 
listic, i.e., a'^{z) - Az^ Q. 

For intermediate times A^^ <C z ^ A^^ photons 
move in the regime of size ~ i/i^T/A , where the pe- 
riodic structure is not resolved and the dielectric cylin- 
ders thereby resemble random scatterers. As usual while 
two optical paths diffusively propagate they find a sig- 
nificant probability to form a common loop and counter- 
propagate along it. Such constructive interference sup- 
presses the normal diffusion coefficient Dd and the weak 
localization correction results, which is order of {vDci)~^ 
with V the photon density of states. Accordingly, the 



linear dispersion is suppressed. That is 

a\z) ^ 
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A-i <z< A-i^ 
> A-i . 



(2) 



z 



We now turn to the detailed analysis of wave dynam- 
ics described by Eq. ^ and derive the rigorous results. 
Below a finite system composed of A/" unit cells will be 
considered with the periodic boundary condition imple- 
mented. In the final results the limit: Af — > -|-oo is set. 

Field-theoretic formalism — First we present the gen- 
eral formalism used throughout this work. Consider 
the propagation on the lattice from 1 to 1' , described 
by the so-called two-point correlator: 3^(1 — = 
N §dvdv'Gf_^^ (r -I- 1, r' -f 1') Gf_^ (r' + 1', r + 1) , notic- 
ing that r , r' are the coordinates in the unit cell. Here 
qR/a jg ^Yi^ retarded/ advanced Green function of Eq. (jT]) 
in the frequency domain. The translational symmetry 
implies that y depends merely on 1 — 1' . 

Following Ref. [l^l the Green functions are written 
down explicitly in the Bloch basis |nk) . In doing so 
the wave dynamics is reduced into the one in a two- 
torus. For the latter we invoke the Wigner transform 
and pass to the phase space T . The coordinate is denoted 
as X = (r, p) (after appropriate rescaling) and satisfies 
Jrj, dX — I . Then, it is a canonical procedure [1, [13, Hi 
to express the spatial Fourier transform of y{\ — l';w) , 
denoted as 3^(q, w) , in terms of the functional integral 
over the supermatrix field Q ~ TAT^^ , where both Q 
and T depend on X ._Moreover, Q_{X) = KC'^QT(X)CK 
and rt(A) = CTT(A:)C'^ with X = (r, -p) . After te- 
dious calculations we find (Throughout this work we are 
interested in the low- lying modes, i.e., q .) 



3^(q,^) 



64A/' 



^ JJdXdX'e''i<''-'''^jDQe 



-F[Q,8i] 



xstr [fc(f + A)(l - T3)g(X)fc(I - A)(I - T3)Q{X')] , (3) 



where = Skj_( 
and the action is 



;.<5k+-k_,q, str is the supertrace, 



F[Q, a]^ dX str TAC^T-^ + —AQ . (4) 



Here all the constant supermatrices: A , ra , fc , K and C , 
as well as the supermatrix fields Q and T are 8x8 matri- 
ces defined on the retarded/advanced, bosonic/fermonic 
and time-reversal sector (for the details of which we refer 
to Ref. i). In Eq. © £a = P • (9r + [zarg, ]) - drV ■ dp 
is the covariant Liouvillian, where the supermatrix a = 
diag(k+, k-)'''' with the superscript standing for the re- 
tarded/advanced sector. Eqs. ^ and (|4|) are the exact 
starting point of the succeeding analysis. 

We remark that the coarse graining T , with the ele- 
ment commensurating with the uncertainty intrinsic to 
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wave dynamics [l^], suffices to regularize the field the- 
ory. (The results below, however, are insensitive to the 
details which, therefore, are not presented.) Indeed, 
the Liouvillian Cq induces a deterministic flow on T : 
X ^ X{t) = (r(t),p(t)). The flow (i) has positive A 
and is ergodic and, (ii) enjoys exponential decay of cor- 
relations and the central limiting theorem for sufficiently 
regular (Holder-continuous) functions. They lead to sin- 
gular hydrodynamic (steady) states (and justify the 
application of the BGS conjecture in the qualitative dis- 
cussions), as we also show below. Then, associated with 
the coarse graining T the Gibb's entropy production re- 
sults which can be shown to be positive-definite [13] , sig- 
naling a stable field theory. 

Notice that properties (i) and (ii) follow directly from 
the choice of the potential V [^. Indeed, the former is 
given by the Donnay-Liverani theorem [14| . To justify 
the latter we introduce the Birkhoff coordinate: {9, ip) , 
where 9 is the position of the incident ray at the dielec- 
tric cylinder boundary {r = rg), and ip £ [0,7r] is the 
angle made between the incident ray and the tangential 
direction of the cylinder boundary. Then, we observe 
that the rotation function A9{(p) defined in [l^, as well 
as its first-order derivative, is discontinuous at (/3 = ^ 
while continuous at each interval. From the Balint-Toth 



theorem 1^ property (ii) immediately follows. 

Hydrodynamic relaxation and Gaspard nonequilibrium 
steady state — Let us start from the ray optics limit: 
A/A — > and present a field-theoretic proof of the de- 
terministic diffusion. The emergence of hydrodynamic 
relaxation from the Liouvillian dynamics is a long stand- 
ing problem ^IG] . It has been exactly solved for the peri- 
odic Lorentz gases as discussed here. To calculate 
Eq. ([3]) we employ the so-called rational parametriza- 
tion: T = 1 + iW . Here WA + AW = and, moreover, 
W{X) = -KCM^^(X)C'^K. Inserting it into the action 
F the perturbative expansion in W shows that is or- 
der of A/ A. Because of A/A in the expansion of 
both F and the prefactor only the quadratic terms are to 
be kept. Upon fixing the C/(l) gauge the succeeding cal- 
culations are formally parallel to disordered systems [^], 
but, there is a crucial difference namely the resolvent: 
{—iu!'^ + £o + P • ig}~^ ■ Recall that, in disordered sys- 
tems, this field propagator is replaced by the heat kernel 
recovering the diffusion. Upon passing to the time do- 
main the resolvent, denoted as 2?g^t , solves 

^s,t = e-*^° - I dTe-(*-")^»p ■ ^gPg,r . (5) 
Jo 

For low-lying modes the functional equation ([5]) gives 



r>g,t = |l - y" dTp(-r) • ig 



(6) 



- dr dr'gg : p(-T)p(-r - r') e 
Jo Jo 



with the second-order iteration. Combined with property 
(i) Eq. ([6]) shows that, given a generic initial distribution, 
for vanishing g it converges to the invariant microcanoni- 
cal measure in the limit t — > -|-oo . A small concentration 
gradient ig drives the system out of equilibrium: 



V, 



= 1 - 



+ 00 



dTp(-T) ■lg + o(g2), (7) 



where the second term is the singular microscopic 
nonequilibrium steady state found by Gaspard [13]. Fur- 
thermore, by inserting Eq. ([6]) into Eq. ((3]) and return- 
ing back to the frequency domain, we find 3^(q, w) = 
2711^/ {—iu! + Dciq^) and, with the help of property (ii), 
rigorously justify the Green-Kubo formula: 



Dr 



- / dX 



+ 00 



dtp{t)-p{0). 



(8) 



Eqs. ^ and ([8]) confirm that the normal diffusion ap- 
pears as the PoUicott-Ruelle resonance as shown ear- 
lier [l^. There, a completely different technique-the 
(■-function theory-is used informing its far-reaching re- 
lation to the field-theoretic approach. In obtaining the 
diffusive two-point correlator the integration over X and 
X' is essential by which the scaling regime is recognized. 

Weak localization — For A/ A ^ 1 wave interference 
accounts for fiuctuations. Keeping the H^-expansion 
up to the fourth order the resolvent is modified to be 
{—iu}^ + Co + 6jC + p ■ ig}^^ ■ Here the perturbed oper- 
ator SC is explicitly expressed as 



6C 



XX' 



dte-"^ \X\V^^\X) 



x{{X\X') - 6{X - X')} 



(9) 



which, in the ray optics limit, is nullified by the fine- 
grained phase space structure because the bases {X\ and 
\X') satisfy {X\X') = 5{X - X') . In wave dynamics 
such point-like structure is smeared by the uncertainty. 
Indeed, a famous theorem of periodic Lorentz gases sug- 
gests that the bases \X) / {X\ must instead be constructed 
as a distribution supported by a small neighborhood of 
X and normalized to unity, with {X\X') defined as the 
overlap between two coarse grained representing points 
X and X' [17]. (Notice that the point-like distribu- 
tion cannot be spanned by these bases.) Furthermore, 
it was rigorously shown [17[ that \X) relaxes to the mi- 
crocanonical measure at some short time scale A"^) 
when the support of \X) uniformly covers T . Upon com- 
bined with Eq. ^ the matrix element (X|Dg_i|X), in 
turn, acquires a universal value while vanishes at shorter 
times t ^ A~^. (For chaotic cavities this fact was noticed 
by Miiller and Altland [l3]-) Importantly, such an inter- 
ference correction conserves the probability (obeying the 
so-called Ward identity) guaranteed by J^dX{{X\X') — 
5(X - X')) = dX'{{X\X') - 5{X - X')} = . 
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Repeat the same procedure of Eqs. ©-([8]) taking into 
account the third order iteration. After tedious calcula- 
tions we arrive at a similar diffusive two-point correlator 
but, crucially, is renormalized into 



1 



1 



(2^ 



1 



(10) 



for A < <C A . The second term fully agrees with ear- 
lier theoretical prediction for periodic Lorentz gases 
and resembles the well-known weak localization correc- 
tion of disordered systems. In the latter systems strong 
localization develops for lower uj . Nevertheless here such 
a scenario breaks down which we now come to study. 

Diffusive-ballistic transport crossover — For lu X the 
zero mode (uniform component) of the T field, denoted as 
To is established recognizing the microcanonical measure 
over T . Therefore, we factorize T ds T ~ T^Ty and 
insert it into the action, where T> are the fluctuating 
component. With T> integrated out Eq. ([3]) gives 



3^(q^0,o.) 



xstr [fc(l + A)(l - T3)Qofc(l - A)(l - t^)Qq] ,(11) 



where Qo 



TqKT^^ and the zero mode action: 



F[go,a] = ^str{i?*[iar3,Qo]' 



-AQo} . Here D* 



ReDo(max{ct', A}) with max{cj,A} playing the role of 
the infrared cutoff. 

To proceed further we invoke the Efetov parametriza- 
tion for Qg 0- Then, the hydrodynamic (i.e., a-) expan- 
sion for Eq. (|lip is performed and kept up to the second 
order. The remaining integral can be exactly done giving 



Kq^O,o.) = 
D{ijj) = D* 



2-KV 



—iuj + I?(a;)q2 
A2 



< A, 



1 



■(1 



2jTrw/A 



(12) 



Strikingly, the diffusive-ballistic transport crossover is 
universal despite that the present unit cell drastically dif- 
fers from a disordered one [l^ . The fundamental struc- 
ture difference at the unit cell level profoundly affects the 
ray optics and thereby the field theory construction. No- 
tice that Del is now renormalized into D* which is not 
reported for periodic disordered systems (jj . These low- 
lying modes are diffusive for sufficiently large frequencies, 
i.e., A oj < X , while ballistic in the opposite limit, i.e., 

< A because of D{u!) w Do(A)A/(-i7rtt') . 

Possible experimental observations — Eqs. ([7]), ([5]), (fTU]) 
and (fT2|) are the main results of this work. To confirm 
the universal diffusive-ballistic crossover might be within 
the reach of the experimental scope of Ref. [l| utilizing 
the optical induction technique [19| . The incident probe 
beam is of 10 ^m full-width at half-maximum. The re- 
fractive index is required to vary over a scale of 100 /im 



which, we expect, might be possible by tuning the op- 
tical interference pattern. The output light intensity 
profile I{x, y; z) then determines the dispersion through 
cr^(z) = J dxdy {x^ + y'^)I{x, y; z) . The latter is directly 
transferred into D{uj) by 



a\z)-a^iO) = J^ 



dcj 1 — I 



Diio), (13) 



giving Eq. ^ with Eq. ((T^) inserted. Finally, we stress 
that the universality of high-frequency transport relies on 
the normal diffusion in the ray optical limit. Thus, the 
predictions are applicable for a large class of periodic di- 
electric structures. In particular, the weak localization is 
stronger for the quasi-one-dimensional transverse plane. 
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